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Abstract 

In this work we propose an action to describe diffusion limited clremical 
reactions belonging to various classes of universality. This action is treated 
through Thompson's approach and can encompass both cases where we have 
segregation as in the A + B ^ reaction, as well the simplest one, namely 
the A + A —> reaction .Our results for long time and long wavelength 
behaviors of the species concentrations and reaction rates agree with exact 
results of Pcliti ioi A + A ^ reaction and rigorous results of Bramson 
and Lebowitz for A + B — > reaction, with equal initial concentrations. The 
different classes of universality are reflected by the obtained upper critical 
dimensions varying continuously from dc = 2 in the first case to dc = 4 in the 
last one. Just at the upper critical dimensions we find universal logarithmic 
corrections to the mean field behavior. 



1 Introduction 

Since the pioneering work of Ovchinnikov and Zeldovich (O Z)[^, followed 
by Toussaint and Wilczek 0, the study of the chemical reactions between 
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two diffusing species which annihilate each other or combine to form an inert 
product when they meet, has becomming a subject of increasing interest [1 
to 8]. As has been recognized since O. Z. work[|T| (see also [2] and [3]), the 
spatial fluctuations in the particle density plays a basic role in the long time 
behavior of their concentrations. Then classical kinetics of reactions rates 
become the proper description of this matter, only above some upper critical 
dimension. 

Early treatments considered situations starting with equal or unequal 
homogeneous concentrations of the two reactant species [1 to 8], but more 
recent studies have also been adressed to the case of inhomogeneous initial 
conditions [9 to 11]. 

A simple example of diffusion limited reaction is given by A+A (inert 
product). A field theory description of this reaction was proposed by Pehti 
[12], performing the renormalization of it at all orders in the perturbation 
theory. 

Krug [13] has proposed a continuous version of a diffusion limited anni- 
hilation (DLA) reaction with a point source. Based on the work of Krug 
[13], one of the present authors [14] wrote na action that when treated by 
the method of the dimensions due to Thompson [15], reproduced the exact 
results of Peliti [12] for the A + A reaction. Thompson's method [15] 
could be thougth as an alternative approach to the renormalization group 
formalism [16]. 

Some important results obtained by Peliti [12] is that the long time be- 
havior of the concentration of species A is given by 



where dc = 2 is the upper critical dimension, namely the dimension above 
which the system behaves classically. 

On the other hand Bramson and Lebowitz[5] performed a rigorous treat- 
ment looking for the long time behavior of the concentrations (e^i) and (e^) 
of two different species in the reaction A + B 0, both for equal and un- 
equal initial concentrations of the two species. In the case of equal initial 
concentrations, namely when ( eA)o = {^b)o, they found that: 




t-i, d<2 
^, d = 2 



t-\ d>2 



(1) 
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being dc — 4: the upper critical dimension in this kind of two species annihi- 
lation reaction. 

A heuristic method as a means to obtain scaling laws for A + A ^ and 
A -\- B ^ reactions was developed by Lindenberg, Shen and Kopelman 
(LSK [17]). LSK work was concerned with the study of the above reactions 
both in euclidian and in fractal geometries and they found that 

(e) ~ t-'^ (3) 

where 



ds 



^ ids 



(4) 



In (4), df and ds are respectively fractal and spectral dimensions, and the 
proper choise of 7 leads to the description of a specific kind of reaction. As 
example, in euclidian geometries, putting dg — df = d, we must choose 7=1 
for the A + B ^ reaction and 7 = for the A + A ^ reaction. 

The aim of this paper is to apply Thompson's method [15] to study the 
long time behavior of some different kinds of diffusion limited reactions. 

C. J. Thompson [15] proposed a simple heuristic method as a means 
to study the critical behavior of a system undergoing second order phase 
transition. He started from a Landau- Ginsburg- Wilson free energy or hamil- 
tonian, and was able to get an explicit relation for the correlation length 
critical exponent (u) as a function of the lattice dimensionality (d). If we 
think that this theory is within the same class of universahty of the Ising 
model, Thompson's work reproduces the exact results for — 2) — 1 and 
= 1) — > 00. 

Thompson's method has been applied to obtain the correlation length 
critical exponent of the Random Field Ising Model by Aharony, Imry and 
Ma [18] and by one of the present authors [19]. His method was also used to 
evaluate the correlation critical exponent of the N- vector Model [20] . Yang- 
Lee Edge Singularity Critical Exponents [21] has been also studied by this 
method. In section 2, we propose an action to describe the A + B ^ 
reaction and we treat it through the Thompson 's approach. This action has 
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been evaluated as an extension of that employed in a previous work [14] in 
order to treat the A +A reaction, and in doing this, an important novel 
ingredient to be considered is segregation. As we will see our treatment is 
able to reproduce some of Bramson and Lebowitz results [5]. In section 3 
we pursue further on this subject, and propose a generalized action which 
reproduces the achievements of section 2 and those of the previous work 
[14] as particular cases. We also will show that this generalization could 
encompass more general situations as some of those described in LSK [17] 
paper. Finally we conclude in section 4 with some comments and prospects. 



In this section we propose to apply Thompson's approach as a means to 
study the two species annihilation reaction (A + S — >^ 0), with equal initial 
concentrations. 

For the case of two species annihilation, at first sight, we would need two 

differential equations to describe its chemical kinetics. However, for equal 
initial densities, due to the simmetry of the problem, the two equations merge 
into a single one. Of course, if the concentrations are equals at starting, they 
remain this way during the passage of the time. 

The starting point of this work is the differential equation: 



where e is the concentration of species(A or B), D is the diffusion constant, 
h{e) stands for an effective input rate and K is the reaction rate. The novelty 
presented by equation (5) as compared to the continuous approximation of 
the A +A — > reaction introduced by Krug [13], is the presence of the 
effective input rate h {t). This term has been introduced in the present work 
as a way to account for the phenomenon of segregation, because we know 
that the mean value (e) decays with the time in such a way that the effective 
input rate h{e) also decreases in time. 

As it is well known, segregation appears in these reactions due to the 
impossibihty of self-annihilation between A (or B) species. 

In order to treat the phenomena described by (5) using Thompson's ap- 
proach, let us define the action: 



2 The A + B 



Reaction. 





(5) 



dt 
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/ 



(6) 



Equation (5) can be obtained from (6), by imposing 5 A — 0, for constant t. 

In an analogous way to the Thompson's reasoning [15] we make the fol- 
lowing heuristic assumptions: 

a- When the integral in (6) is taken over the cube l"^ in d dimensions, the 
absolute values of the four terms separately in (6) are all of order unity. 

b- K{1) is finite in the limit / — > cxo. 

Using assumption (a), that is a scaling assumption, in the first term of (6) 
we have: 



/ d'^r 

Jld 



f-'{e') ~ 1, 



so that the mean squared value of e ((e^)) behaves as 

/ 2\ i2-d 

For the fourth term in (6), we have 



f d'^r 

Jld 



2 dt 



where we have considered that 

e {t) = eo cxp {~Qt) 
Putting (8) into (9), we obtain: 



(7) 



(8) 



(9) 



(10) 



(11) 



We observe that (11) points out to the signature of the brownian character 
of these diffusion limited reactions. 

Applying hypothesis (a) to the second term of (6) leads to: 



Id 



d'^r{h{e)e] 



2 id 



{h){eyi 



(12) 
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Supposing (/i) ~ 1, we have for the mean value of e ((e )) the behavior: 

(e) ~ H (13) 
Using (11) into (13), we obtain: 

(e) ~ (;')~^ ~r-i (14) 

This result (Eq. (14)) agrees with that obtained by Bramson and Lebowitz 
[5] which performed rigorous calculations for the A + B reactions. We 
observe that we have treated here only the case of equal initial concentrations, 
namely (e^)o = (es)o- The case of unequal initial concentrations will be 
treated elsewhere. 

Now let us use assumption (a) in the third term of (6), 

/^A(life») 

In order to pursue further let us make a plausible hypothesis: 

(6^) ^ {e'){e) (16) 

This kind of decouphng could be justified taking in account that the reaction 
term which appears in Eq. (5) is essentially of a bilinear form ( a two points 
correlation). Beside this we observe that this decoupling has been worked 
well in a previous calculation [14], where we reproduced the exact results of 
Peliti [12] for the A + A^O reaction. 

Using (8), (13) and (16) in (15), we obtain: 

(K) ^ 1-2-^ (ff'' (17) 

By considering (11), we can write also: 

{K) ~ T^-i (18) 

Equation (17) displays d — 4, ets the upper critical dimension for the model. 

We see from ( 17) that ( K) diverges as / — * oo, for d > 4. Then, in order 
to satisfy assumption (b) we adopt ( K) = 1, for d > 4. So the mean field 
description of this model is correct for d > 4. 



3\ jd 



(15) 
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Now, at the upper critical dimension {d = 4), let us make some additional 
considerations in order to obtain the logarithmic corrections to the reactant 
concentration e. We have for d — 4:( see (8), (13) and (16)): 

(e^) - (e2)(e) for d = A. (19) 

If wc use assumption (a) for the third term of (6) in a somewhat modified 
context, and after having made inside the integral the following substitution: 

~ r-^ (20) 

we can write 

j^^ d'^r (^^e^) ~ {\^^~^) ^^^^ ^ ~ ^ 

where / is measured in units of space lattice, being 1 the lower cutoff. 
Then we have for d = 4, 

{K)^Mir\ (22) 
At this point, let us consider the following differential equation: 

This equation can be obtained by considering the equality of the first and 
third terms of (6)( given by assumption a), where we suppose that we can 
substitute equality between integrals by equality between integrands, after 
replacing K by {K). 

We can solve ( 23) at the upper critical dimension {d — A),hy performing 
the integration between 1 and /, and using (22). We get: 

(,) . . Mi), f„ rf = 4. (24) 

t T" 

As we can see from (22) just at the upper critical dimension for the A + B ^ 
reaction {dc = 4), we have obtained a logarithmic correction to the mean 
field description. The same behavior was obtained before for the A + A — > 
reaction at dc = 2. (see [14]). 
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3 Unified Treatment oiA+A^O and A+B ^ 
reactions. 

The relative success of Thompson's approach [15] in treating both A+A — > 
and A + B ^ reactions encourage us to look for an extended action which 
could encompass the two mentioned reactions as limiting cases. So let us 
write the following extended action. 



A^= / d'^r 



2 ' ' ' ' 3 2 dt 



(25) 



where cr > 0. 

In an analogous ways that we have worked before, we can apply Thomp- 
son's method [15] to the action giving by (25), obtaining: 

(e),~r'^/(^+^), (26) 

and 

{K),^[f)^''\ (27) 

Relation (27) imphes that, for the model defined in (25), the upper critical 
dimension is given by 

4 {(t) = 2(7 + 2. (28) 
So alternatively we can write {K)^ as 

{K),^(l^Y^'' r^A-\ (29) 

We observe that \i a — dc — 2 and we recover the A-\- A ^ behavior 
[12,14], and putting o" = 1 we get dc = 4, reproducing the results of the 
previous section for the A + B ^ reactions. 

It is interesting to verify that the action given by (25) can reproduce the 
results of the LSK work [17], in the case of euclidian geometries. First we 
can write: 

(e).-r-'='/(2-+2) ^^-d/dc, (30) 

Making the requirement that the exponent of right side of (30) to be equal 
to —a given by (4), and by setting ds = df = d, we have: 
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Therefore putting 7 = a = 0, and 7 = a = 1, we get the long time behavior 
of the A + A ^ and A + B ^ reactions respectively. The case of 
A + B — > reactions where correlations are important (see [17]) can be 
described for some intermediate values of the parameters 7 (or a) . 

Finally is worth to point out that at the upper critical dimension of the 
model described by the action (25) we have again the logarithmic corrections 
to the mean field behavior. This result can be obtained in the same way we 
have worked out before in section 2. Doing this, we obtain: 

(e),.!!i».!!iM. for 4 = 2. + 2. (32) 

T 

Before concluding this section we would like to write an extendend version 
of the renormalized reaction rate in the same spirit of Peliti's work [12]. For 
long time and long wavelengths, we can write the scaling relation: 

(33) 

Putting o" = 0, we recover the renormalized reaction rate of Peliti [12] ( 
equation (10) of his paper). 




4 Concluding remarks. 

In this paper we have proposed a model lagrangian ( action) as a means to 
describe a large class of diffusion limited reactions. The long time behavior 
of the concentration e and of the reaction rate K were evaluated by using 
Thompson's approach [15], which can be thought as being a simple alternative 
way to the renormalization group method [16]. 

In two particular cases of this extended action, we recover the exact re- 
sults of Pehti [12] ior A + A — >■ and the rigorous bounds of Bramson 
and Lebowitz[5] for the A + 5 — > reactions. Besides its simplicity, one of 
the advantages of the present method is that it can encompass diffusion lim- 
ited reactions of different classes of universality within the same formalism. 
Moreover, our calculations display logarithmic behavior of the species con- 
centrations and reactions rates at their respectives upper critical dimensions 
for any class of universality. 
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An other advantage of the present method is that it can be extended to 
study, for example, the scalar field theories and for dimension d, 

in such a way that we can obtain a certain critical dimension dc for a given 
value of n in theory in a closed form. But now, dc must be interpreted 
as a dimension where the theory becomes renormalizable displaying log- 
arithmic dependence of the coupling constant in the energy scale. This line 
of reasoning will be subject of a forthcomming paper. 
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